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Abstract 

Surface catalytic processes produce, under certain conditions, small 
clusters of adsorbed atoms or groups, called islands which, after they 
have been formed, move as individual entities. Here we consider the 
catalytic reduction of NO with hydrogen on platinum, (i) Using video 
field ion microscopy, we observe the dynamic motion of small hydroxyl 
islands on the Pt(OOl) plane; despite changes in their morphology, the 
islands dimensions are confined to values corresponding to 10 to 30 Pt 
atoms suggesting cooperative effects to be in operation, (ii) We con- 
struct an automaton (or lattice Monte-Carlo) model on the basis of a 
set of elementary processes governing the microscopic dynamics. The 
agreement between the simulation results and the experimental obser- 
vations suggests a possible mechanism for the formation and dynamics 
of hydroxyl islands. 

Keywords: Surface catalytic processes, Adsorbate islands. Field ion 
microscopy, Cellular automata. 
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Modem surface science aims at elucidating the microscopic mechanisms of 
dynamic reaction phenomena Scanning tunnehng microscopy (STM) has 
largely contributed to the development of the field by providing information 
at the atomic scale. For example, a recent study of the oxygen dissociation 
on the plane (111) of a platinum surface [@] has made visible the formation of 
adsorbate clusters called islands. Field ion microscopy (FIM) is an alterna- 
tive method where samples are given in the form of small three-dimensional 
tips exposing a number of crystallographically different planes. The capabil- 
ity of the method to image dynamic reactive phenomena at nanoscale has 
been explored 0, ^, |^. For example, reaction-diffusion fronts moving across 
the tip surface have been made visible thus providing information on the 
communication behavior of different planes. As compared to STM, the time 
resolution of the FIM method is only limited by the frequency of the video 
system (20 msec. here). 

While the formation and dynamical behavior of islands appear as a se- 
quence of processes which, in general, are difficult to explain on the basis of 
a phenomenological description, this type of system is well suited for mod- 
eling with Cellular Automata 0, which are constructed with simple rules, 
and may thereby suggest basic mechanisms for complex reactive phenomena. 
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Here we report on FIM experiments showing the formation of small clusters 
of hydroxyl species from a co-adsorbed 0^^ /^ad layer on the (001) plane of 
a Pt tip surface during the reaction of nitric oxide (NO) with hydrogen gas, 
and, in order to gain insight into the basic mechanisms governing the fairly 
complex processes producing islands, we construct an automaton (or lattice 
Monte-Carlo model) with simple rules which offers a microscopic approach 
applicable to the experimental system studied by FIM. 

The experimental studies were performed with a [001] -oriented Pt tip 
(metal purity 99.99 %) prepared by electrochemical methods. The sample 
was mounted in an all-metal ultra-high vacuum (UHV) field ion microscope 
(base pressure 10^*^ Pa) and cleaned using standard methods |^, 0. Gases 
were used in either commercially available purity (NO: 99.5 %) or after fur- 
ther purification using adsorption methods (H2 and Ne: better than 99.999 % 
each). Ne-field ion imaging prior to reaction experiments indicated a clean 
defect-free Pt tip. As described in the nearly hemispherical shape of the 
tip specimen is transformed into a pyramid during adsorption of NO gas or 
reaction with NO/H2 gas mixtures, respectively. An essential feature of the 
pyramidal form is the appearance of large (111) oriented slopes and a (001) 
truncated top considerably increased in size as compared to the original state. 
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Figure 1 shows a respective Ne-field ion image at 57 K obtained after the 
reaction studies. Only a few net planes are observed in the region between 
the (001) pole and the peripheral (111) planes. Subsequent to NO/H2 re- 
action studies, a chemical surface analysis was performed using atom-probe 
techniques. The respective mass spectra P| demonstrated the Pt tip sur- 
face to be oxidized, i.e. Oa^-covered. While in previous research, oscillating 
reaction phenomena with water (and nitrogen) product formation were ex- 
amined in detail |^ , we focus here on the occurrence of OHq^ islands from 
an Oad /Had co-adsorbed layer, with particular emphasis on the dynamics of 
island formation and displacement on the Pt (001) pole. 

The studies were performed under truly in-situ conditions, i.e. video field 
ion imaging during the ongoing surface reaction. Figure 2 shows a typi- 
cal sequence of images as obtained while exposing an initially clean, nearly 
hemispherical Pt tip specimen to reactant pressures of p//2=4xlO^^ Pa and 
PAfo=3xlO~'^ Pa at 500 K and 8.7 V/nm. After reaction-induced shape 
transformation into a pyramidal morphology small islands (appearing as dark 
patches on a bright background in Fig.2) with an equivalent size of 10 to 30 
Pt surface atoms become discernible on the (001) pole. Islands exclusively 
form at the layer edge and move into the (001) terrace region. Up to four 
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Figure 1: Field ion micrograph (image gas: H2) of a Ft tip demonstrating 
a truncated pyramidal form with a [001] oriented top; the morphology is 
the result of a NO-induced reconstruction of an originally hemispherical Pt 
specimen; the image was obtained after field evaporation of a few surface 
layers subsequent to reaction with NO/H2 gas. 
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islands can be seen at the same time. Mean lifetimes of several minutes are 
observed for individual islands before their annihilation at the layer ledge. 

One of the interesting features in the genesis of small islands is their crit- 
ical size. Accordingly, clusters extending over less than ~ 10 atoms of the 
Pt(OOl) plane are not observed. Conversely, spreading over more than ~ 
30 Pt atoms does not occur either. Moving clusters yet suffer some changes 
in their morphology. Their simultaneous presence on the (001) plane may 
cause collisions, however, no mergence is seen. Under the experimental con- 
ditions applied, both NO and H2 undergo dissociation on the Pt tip sur- 
face. The dissociation probability is strongly dependent on the local surface 
crystallography. Accordingly, starting with a clean surface, the Pt (001) 
plane is readily covered by O^^ (nitrogen atoms formed concomitantly dur- 
ing NO decomposition undergo recombinative thermal desorption at the re- 



action temperature, T=500 K [[12], [T^, |Tj]); while the rough surface planes 



accommodate both Oad and Had ■ In fact, hydrogen activation at the (001) 
layer edge must be regarded a critical step preceding Had injection into the 
Oad covered (001) plane. Following the resonant field ionisation mechanism 



proposed by Kreuzer and Wang |T5[, an oxygen- covered Pt surface is ex- 



pected to be imaged with considerable brightness when using molecular NO 




Figure 2: Typical sequence of video field ion images demonstrating OHa^ 
island formation during the catalytic reduction of NO by H2- Experimental 
conditions: T = 500 K, p^^^ = 4 x lO'^ Pa, pjvo = 3 x 10-^ Pa, F= 8.5 
V/nm. Micrographs are treated by software allowing to optimize contrast and 
brightness; circles are understood as a guide to the eye to approximate the 
location of the (001) plane. Snapshots taken over a total time interval of order 
of one sec. demonstrate islands displacement, collision and separation; the 
arrow points to an island that will dissolve at the layer edge while proceeding 
to (f). 
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gas. On the other hand, hydroxy! species formed by reaction of Oad with Had , 
are thought to be imaged with less brightness than 0^^ , mainly because of 
the lower NO ionisation probability at slightly larger critical distances from 
the surface. The above arguments lead to the conclusion that OHq^ species 
are formed during the NO/H2 reaction. They appear as small clusters with 
low FIM brightness in the "defect" region of the (001) layer edge and move 
collectively into the flat terrace region. One of the attractive features of 
the underlying reaction mechanism is the displacement of hydrogen islands 
on top of Oad rather than exchanging Oa^ for OUad ■ Alternative explana- 
tions according to which field-adsorbed water clusters are associated with 
the occurrence of FIM islands, can be rejected on the basis of results ob- 
tained for the O2/H2 reaction on Pt tips by field ion appearance potential 
spectroscopy 

So far the microscopic mechanisms governing the dynamics of the exper- 
imentally observed phenomena described above have not been given a satis- 
factory explanation. A possible approach to the microscopic aspects of these 
surface reactive phenomena is provided by cellular automaton approaches 
or lattice Monte-Carlo modeling |T^. Here we construct a cellular automaton 



[O] on a two-dimensional square lattice whose nodes can be occupied by vir- 



tual particles with exclusion principle (no more than one particle per node). 
The particles are subject to probabilistic motion with displacements along 
the four lattice directions: at each time step each particle can hop to any of 
its four nearest neighboring sites. The exclusion principle precludes simul- 
taneous displacements of all the particles in the automaton universe in one 
time step. Therefore a sequential updating is implemented: the propagation 
phase consists of a sequence of successive displacements (on the average one 
per particle) and the effective time step is defined as the sum of N automaton 
time steps, where N is the number of particles in the automaton universe. 
Each particle can move to one of its neighboring sites, unless its destination 
site is occupied, in which case the particle does not move. The particles are 
subject to non-local two-body interactions, which are short-range attractive 
and long-range repulsive, and the amplitude and the range of the attractive 
well and of the repulsive wall can be tuned parametrically. 

The justification of the automaton for the modehng of the surface phe- 
nomena observed in the experiments showing adsorbate islands, is as follows: 
(i) the processes take place in two-dimensional space, the Pt surface, and 
the (001) top plane has square symmetry (bulk-truncate form); the automa- 
ton reproduces these conditions as well as the size of the top plane (30 sites 
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across); (ii) the restriction of a single hydroxyl group per site is accounted 
for by the exclusion principle; (iii) an hydroxyl group can be identified as 
an automaton particle; (iv) the automaton has intrinsic fluctuations; (v) the 
particle dynamics depends on the interactions between OH groups (see cap- 



tion of Fig. 3) which are short-range attractive and long-range repulsive [|T£ 

In the "vacuum" a particle has equal probability to move to any of its four 
neighboring sites; in the presence of other particles, this probability is biased 
according to the occupation of the sites within some neighborhood around the 
particle. The neighborhood is defined as the lattice domain surrounding an 
occupied site and containing all sites within interaction range. The bias is 
such that the sum of the probabilities at each site is conserved (= 1) so that 
a positive (negative) interaction along one channel induces a repulsive (at- 
tractive) interaction along the other channels, and the propagation direction 
is drawn according to the probabilities evaluated, for the different channels, 
from the occupation of the lattice sites in the neighborhood. 

The updating of the asynchronous automaton at each effective time step 
consists of a sequential series of operations, (i) An occupied site is chosen 
randomly; the probability that the particle be displaced to one of its four 
neighboring sites is initially set to .25 for each channel, (ii) The sites located 
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within the neighborhood of the chosen particle are tested for the presence of 
particles, and the displacement probabilities are adjusted according to the 
position and distance of the particles within interaction range, (iii) In order 
to preserve the stochastic nature of the procedure, a random number is drawn 
between and 1 and compared to the cumulated values of the probabilities 
of the successive channels; the channel for which the cumulated value exceeds 
the random number, is selected, (iv) The particle is displaced by one lattice 
unit along the selected channel, unless the destination site is occupied, (v) If 
so, the particle does not move, and the propagation probability is transfered 
to the target particle which memorizes one "coUisional unit" in the corre- 
sponding channel for later updating; the initial probability for displacement 
(.25) is then biased accordingly. This step in the algorithm acts as a process 
with momentum transfer memory. 

The boundary conditions of the automaton universe should be such that 
they correspond to the physical situation where hydrogen invades the Pt(OOl) 
surface from the adjacent planes, and particular properties must be assigned 
to the border nodes: a particle located on a site within a given distance 
of a border node is subject to a repulsive force which favors the particle 
displacement along the channels pointing inwards the system. This affects 
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Figure 3: Islands formation and dynamics in automaton with circular open 
boundary (diameter: 30 nodes). Interaction range and amplitudes: 9 nodes, 
with attraction on first 2 nodes and next 3 nodes, zero-interaction over next 2 
nodes, and repulsion over last 2 nodes; relative values of the amplitudes: —6 : 
—3 : : +2 respectively. Repulsion from boundary nodes: interaction range 
9 nodes, with relative amplitude corresponding to 5% bias along the inward 
direction. Burst injection rate: 30 time steps every 100 time steps. The 
time sequence goes from left to right and from top row to bottom row: time 
t = 1 (upper left panel), 600,1500,3300,5100, and 5700 (lower right panel) 
{t in automaton time step units). An animated version of the simulation is 
available on the web site |http : / /poseidon . ulb . ac . be . 

any particle residing on a site located at a distance equal or smaller than the 
border interaction range. Circular boundary conditions are imposed (i.e. a 
set of nodes with the boundary 

repulsive effect is chosen to approach a circular geometry on the square lat- 
tice), and particles are injected in the form of bursts over a limited region of 
the boundary as follows. The circular boundary is ascribed a given arbitrary 
angular reference location where from an angle a is chosen randomly. The 
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set of boundary sites in the arc spanned between a and a randomly chosen 
Aq; define the region where particles are injected randomly during a time At; 
the time intervals between Ar's are periodically distributed. The injection 
procedure is repeated for the duration of the simulation. 

In Fig. 3 we present a sequence of snapshots from the model simulation 
starting with a burst injection and leading to the formation of islands and 
their subsequent evolution. Comparison with the experimental micrographs 
of Fig. 2 shows that the clusters obtained in the automaton simulation ex- 
hibit the essential features of the observed islands: their size (~ 30 sites) 
and number (3 to 4) have the correct order of magnitude, they remain "au- 
tonomous" (when in contact they do not merge) , and their dynamics reflects 
the erratic motion of the islands. The agreement between the simulation 
results and the experimental observations suggests that the elementary pro- 
cesses on the basis of which the automaton is constructed, provide a plausible 
mechanism for island formation and dynamics as observed in FIM experi- 
ments whose results were awaiting satisfactory explanation. The basic mech- 
anism responsible for the islands formation is governed the balance between 
attractive and repulsive forces between OHa,^ groups. The relative values of 
the interaction forces given in the caption of Fig. 3 give indication as to what 
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a plausible interaction potential could be between OHa^ groups; furthermore 
we found the dynamics to be quite sensitive to the value of the parameters 
(for further explorations of surface aggregation phenomena by automaton 
simulations, see |^), and the values used here optimize the agreement with 
the experimental observations. The other crucial feature is the role played 
by fluctuations as intrinsically contained in the automaton dynamics. So we 
conjecture that island formation and displacements are governed primarily 
(i) by the balance between lateral attractive and repulsive forces between 
OUad groups, and (ii) by dynamical fluctuations with memory effects caused 
by frustrated momentum transfer. 
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